Minimal qudit code for a qubit in the phase damping channel 
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Using the stabilizer formalism we construct the minimal code into a D-dimensional Hilbert space 
(qudit) to protect a qubit against phase damping. The effectiveness of this code is then studied by 
means of input-output fidelity. 
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I. INTRODUCTION 

Quantum Error Correction (QEC) theory [l|, d, H, Q 
concerns the possibihty to protect against environmental 
noise when storing or transmitting quantum information. 
This possibility relies, likewise in the classical error cor- 
rection theory, in redundancy. This implies embedding a 
quantum information unit (a qubit) belonging to a given 
Hilbert space {'H2) into a larger one {Hd with D > 2). 
The latter is usually chosen as n times the tensor prod- 
uct of the former, so that D = 2" and Hd = nf". We 
generally refer to this kind of encoding 



(1) 



as block- encoding of the qubit. However, it is obvious 
that whatever Hd (with D > 2), extra state space is 
available and could potentially be exploited without any 
restriction. As a consequence, the alternative possibil- 
ity is to embed a qubit into a _D-dimcnsional quantum 
system, i.e., a qudit, where D can be made even infinite 
(in the limit where the qudit becomes a quantum oscil- 
lator, i.e., a bosonic mode [1,0). We refer to this kind 
of encoding 



Hd H 



(8n 



(2) 



as to qudit- encoding of the qubit. In a standard QEC 
framework, such encodings are coupled to suitable decod- 
ing stages, where a recovery operation (e.g., syndrome ex- 
traction and error correction) restores the original quan- 
tum information by removing the (correctable) errors in- 
duced by the noisy action of the environment. 

Note that the encodings of Eqs. ^ and ^ are not 
equivalent. A first simple reason relies in the available 
dimensions for qudit-encoding, which are not necessarily 
restricted to powers of 2. Moreover, at a more fundamen- 
tal level, the errors affecting the two storing systems are 
different, that is, they form two different algebras. For 
a block of qubits errors are given by combinations of bit 
and phase flips, which are representable in terms of prod- 
ucts of Pauli matrices. A single qudit instead is affected 
by amplitude and phase shifts (implying that, asymp- 
totically, a single bosonic mode is affected by diffusion 
in position and momentum Q), which are represented 



by the unitarily generalized Pauli matrices for a single 
qudit. 

Pioneering advances in the QEC with higher- 
dimensional spin systems 0, Q and bosonic modes ^ 
were achieved during the nineties. More recently, Ref. 6| 
introduced novel kinds of codes for qudits, known as shift- 
resistant (SR) quantum codes. In its simplest formula- 
tion, a SR code corresponds to embedding a logical qubit 
into a larger qudit, followed by a recovery stage which 
restores the quantum information from a bounded set of 
quantum errors (i.e., amplitude and phase shifts whose 
weight is less than some critical value). In particular, 
Ref. showed that a qudit of dimension D = 18 rep- 
resents the smallest quantum system able to protect a 
logical qubit from a single quantum error, where the cor- 
responding 5-qubit block code [[5,1,3]] of Ref. needs 
a Hilbert space of dimension £> = 2^ > 18. Let us un- 
derline that both of these codes are stabilizer codes Q 
and are perfect, roughly meaning that they need minimal 
quantum resources for their task [To| . 

The latter peculiarity is very important since the pri- 
mary issue for having experimentally feasible QEC codes 
consists in simplifying their complexity. In fact, the im- 
portance of using minimal resource codes relies on our 
current difficulty in performing high fidelity operations 
on a small number of qubits [111 ]. However, apart from 
the optimality of the above perfect codes (which are de- 
signed to defeat general quantum errors), it is still an 
open problem to find the most efhcient quantum codes 
which enable QEC within .specific error models. In fact, 
if the dominant decoherence process in a physical system 
is of a specific nature and well known, one can look for 
a corresponding quantum error correction scheme whose 
quantum complexity is as small as possible. Such a prob- 
lem has been raised, for the first time, in Ref. for pro- 
tecting logical qubits against dcphasing. Later, Ref. [13] 
proposed an optimal code embedding a qubit in a block 
of bosonic modes able to protect against the effect of 
amplitude damping. 

To date, nobody has analyzed the same problem for 
qudit-encoding, i.e., nobody has considered the engineer- 
ing of a minimal singie-qudit code able to protect a log- 
ical qubit against a specific kind of decoherence. Only 
Ref. [ij] pointed out that qudit-encoding is not offec- 
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tive by itself when specific error models are taken into 
account. That is, without a suitable error correction (re- 
covery) operation, the extra space cannot be exploited 
to protect against errors. In this paper, we consider the 
qudit-encoding in a QEC framework (i.e., with a suitable 
recovery stage) and we design the minimal codes which 
are able to protect a logical qubit against a single class 
of errors, such as amplitude or phase shifts. Note that 
we are here considering minimal codes which are quasi- 
classical. In fact, even if they encode quantum informa- 
tion (one logical qubit), the environmental error-model 
here is classical, in the sense that the correctable errors 
occur only in a preferred basis. The two complementary 
bases of a single qudit are perfectly symmetric, being 
connected by a discrete Fourier transformation. There- 
fore, by fixing the unperturbed basis (pointer basis) to 
be the computational one, we can always define as phase- 
damping the damping that affects the complementary ba- 
sis. We will show the robustness of a minimal qudit code 
in preserving the encoded quantum information against 
this kind of error. 

The layout of the paper is the following. In Section |TT] 
we present the code's construction and its performance 
against shift-errors. Section IIIII is devoted to the phase 
damping channel. Section IIVI shows the performance of 
the code against phase damping in terms of input-output 
fidelity. Finally, Section |V] is the conclusion. 



II. THE CODE 
A. Qudits 

Let us consider a qudit, i.e., a ZJ-dimensional spin- 
system. In its Hilbert space Tio we choose a computa- 
tional basis {\j)} labeled by modular integers j G Z^i := 
{0, . . . , Z? — 1}. An arbitrary unitary transformation 
Ti.u —> Tiu can be expanded in terms of generalized 
Pauli operators 

X^Z^ , a,6G Zi) , (3) 

which are defined by 

X \j) - 1.7 ® 1) , Z \j) = \j) , (4) 

where ji © j2 := ji -I- j2(nrod£') and 

uj cxp(i27r/i:>) . (5) 

Such unitary operators satisfy the anticommutation re- 
lation 



ZX = LoXZ , 
and their eigenstates are connected by 

_ D-l 
i=0 



(6) 



(7) 



where X\i) - 
with entries 



and H is the D x D Fourier matrix 



Hi 



i,j e 



(8) 



Accordingly, a general quantum error acting on the qu- 
dit can be decomposed in the error basis of Eq. ^ . Its 
elements, i.e., the generalized Pauli operators, represent 
the basic quantum errors which a quantum correcting 
code must correct. According to Eq. ([4]) these are dis- 
tinguished as amplitude shifts X° and complementary 
phase shifts Z^. Multiplying by suitable phase factors 
the elements of Eq. ([3]) , one defines the qudit Pauli group 
and, consequently, extends the stabilizer formalism [3|,|6|. 
These two kinds of errors, if considered separately, rep- 
resent an abelian group and their correction can thus be 
performed through quasi-classical codes. 



B. Single errors 

It is natural to ask what is the smallest D-level system 
which protects an encoded qubit from a single amplitude 
shift X^^ . Let us first consider an example with D = 6, 
so that Lu = exp (z7r/3). A logical qubit can be encoded 
in the two codewords stabilized by the generator Z^, i.e.. 



where 



|0) := |0) , |1) := |3) , 



Z^ |0) = |0) , Z2 |3) = |3) 



(9) 



(10) 



In such a case the measurement of the stabilizer preserves 
every coherent superposition \ip) = a |0) + /3 |3), i.e., 

Z'\^) = \^), (11) 

while it detects single X errors, i.e., 

Z^X^^lifi) =LU^'^X^^\ip) . (12) 

Alternately, one must consider the complementary gen- 
erator X^ for correcting single Z^^ errors, i.e., one must 
encode the qubit into the codewords 



where 



1+) := |0) , I-) |3) , 



X^IO) = |0) , X^\3) = |3) 



(13) 



(14) 



According to Eq. ([7]), one can express these codewords in 
the computational basis as 



j=o 



V6 
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C. Multiple errors 

From the previous example we argue that in order to 
correct k shifts we need a qudit with 



D = Ak 



(16) 



levels. This can be understood by means of the "clock" 
picture of Fig. [U and is easily proven in the following. 
Consider, for example, the case of k amplitude shifts 




FIG. 1: Pictorial view of errors' effect on code states. The 
shadowed area separates the correctable error spaces associ- 
ated with the two codewords. 

X^'^ , but the reasoning is perfectly symmetric for the 
complementary errors. In order to determine two possi- 
ble codewords, we must consider the eigenvalue equation 



(17) 



where j =^ G Cz Ijd and D must be determined. 

In Eq. (fT7|) . the state \j) is stabilized if and only if oj^'-^ = 
1, and this happens in either the trivial case j = or in 
the case 



(18) 



Note that, since j < D, we must necessarily set G > 2 
for the weight of the Z-generator. The resulting code 
will be able to correct all errors . . .X^*^ 

if and only if the stabilizer generator Z'~^ will commute 
with X'^, where d := 2k + 1 defines the distance of the 
code. Since 



Z'^X'^ = uj'^'^X'^Z'^ 
this will happen if and only if 

dG = D , 



(19) 



(20) 



i.e., if the weight of the generator corresponds to the 
the ratio between the dimension of the qudit and the 



code's distance. By comparing Eqs. ((T8)) and ([20]) we 
must conclude that 



(21) 



j = d = 2k + l 



Then, by making the "minimal" choice G = 2 in Eq. (|18p. 
we get the minimal dimension D of Eq. p6|) . In con- 
clusion, k amplitude-shifts are corrected by enconding a 
qubit into the codewords 



|0) :=|0) , |l):=|2fc + l) 



(22) 



of a {4k + 2)-dimensional qudit. These codewords are 
stabilized by the operator Z^ and are connected by the 
logical flip gate X :— X'^^'^^. Analogously, k phase-shifts 
are corrected by means of the codewords 



+ : = = 



1 



4fc+l 



V4F 



(23) 



3=0 



4fe+l 



|2fc + l) 



V4F 



which are stabilized by X^ and are connected by the 
logical phase gate Z := Both these codes are 

perfect, since the correctable error spaces (of dimension 
d) associated to their codewords just barely fit in the 
qudit space (of dimension D = 2d). We may refer to 
these codes as to minimal amplitude code {|0) , |2fc + 1)} 

and minimal phase code {|0), |2fc -t- 1)}, respectively. It 
is clear that they are equivalent up to a (discrete) Fourier 
transformation. 

Note that, correspondingly, the minimal qubit block- 
code which is able to correct k phase (or amplitude) error 
flips works via majority voting and, therefore, needs a 
block of 2 A: -|- 1 qubits. This is equivalent to considering 
a Hilbert space of dimension D = 2^'^+^ which is ex- 
ponential rather than polynomial in k. This means that 
the qubit code is exponentially more demanding than the 
corresponding shift-resistant qudit code at given weight 
k. 



D. Syndrome extraction and error recovery 

1. Amplitude errors 

Consider an arbitrary coherent superposition 
of orthogonal codewords of the amplitude code 
{|0),|2fc + l)}, i.e.. 



1^(0)) =a|0)-f /3|2fc + l) 



(25) 



Suppose that an amplitude shift-error X'*, with syndrome 
—k < s < k, occurs on this superposition. Then, the 
logical state of Eq. ((25|) becomes 

\ip{s)):^X''\^(0))=a\0®s)+P\2k + l®s) . (26) 
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According to Sec. Ill C[ such an error is detected by mea- 
suring the complementary generator . In fact, this 
measurement gives 

Z^\^{s))=u?^\^{s)) , (27) 

i.e., the syndrome is unambiguously extracted via the 
eigenvalue uP'^ of (non-degeneracy of the code), while 
the corrupted state is preserved in the process. In order 
to realize this kind of quantum non-demolition measure- 
ment, wc must append an ancillary system to the signal, 
let the joint system evolve according to a suitable uni- 
tary interaction, and finally measure the ancilla. Since 
we must distinguish Ik + 1 orthogonal errors (fc positive 
shifts, k negative shifts, and the no-shift), we need an 
ancillary system having at least 2fc-|- 1 orthogonal states, 
that we label by |^)^ with I G Z, \l\ < k. 

In detail this correction process goes as follows. Let us 
introduce the 2k + 1 projectors 

P{s) |o e s) (0 e s| + \2k + 1 e s) (2fc + 1 e s| , (28) 

and construct the following unitary operation (general- 
ized CNOT) 

k 

N:=Y^ P{s)X% . (29) 

s— — k 

It is then easy to check that N realizes the syndrome 
extraction. In fact, its effect on the joint system signal- 
plus-ancilla corresponds to leave the corrupted state un- 
changed while shifting the ancilla by a quantity equal to 
the syndrome, i.e., 

iV(|(^(s'))®|0)J = |^(s')>®k% . (30) 

At this point, the measurement of the ancilla provides 
the syndrome s' and one restores the original signal state 
by applying the corresponding inverse operator X^" to 

l^(«')>- 

It is known that the last procedure, i.e., the error cor- 
rection stage, can be also implemented in a unitary man- 
ner. In fact, we can define the correction operator 

fc 

^ , (31) 

s— — k 

and applying it to to the final state of Eq. (|5D|) . In such 
a way we get 

c(|^(.'))®k%) = l^(o))®k% , 

thus recovering the initial encoded state of Eq. ((25|) . Note 
that the two unitary operators of Eqs. and PT|) can 
be compacted together in a unique recovery operator 

fc 

R:=CN= ^ X-^P{s)(^\r)^(r\X^^. (32) 

r,s— — k 



The above derivation simply shows how the recovery pro- 
cedure works properly when a logical state |'p(0)) is af- 
fected by amplitude error-shifts which are correctable, 
i.e., which fall within the distance of the code. More 
generally, we may ask how the recovery works when such 
errors are not necessarily correctable. This is a question 
that must be answered if we want to test these codes in 
a quantum communication scenario where the decoher- 
ence of a channel can be very strong. To this purpose we 
must first derive the effect of recovery on a completely 
arbitrary state of the qudit. 

Thus, let us consider an arbitrary state 

D-l 

of a qudit with dimension D = 4fc + 2. The joint action 
of the recovery operator reads 

i?(p®|0)^(0|)i?t 

fe 

= X-^P{s)pP\s')X^' ®\s)j,{s'\ . (34) 

s.s' — — fc 

If we now trace out the ancilla, we get the recovery map 
acting on the qudit state 

fe 

£n{p)=Y.X-^P{s)pP\s)X^. (35) 

s— — fc 

By virtue of Eqs. and it is equal to 

Er{p) = $(0, 0) |0) (0| + $(0, 2k + 1) |0) {2k + 1| 
+$(2fc + l,0)|2A:+l)(0| 
+$(2fc+l,2A: + l)|2fc + l)(2fc + l| , (36) 

where 

fc 

^{x,y) := ^ Px®s,y®s ■ (37) 

s— — fc' 



2. Phase errors 

For correcting phase errors the procedure is perfectly 
analogous to the previous one. It is sufficient to exchange 
the role of X and Z and ^counWor the rotated code- 
words of the phase code {|0), |2fc -f 1)}. Thus, the action 
of recovery on an arbitrary state p of the system is now 
described by the map 

Er{p) = $(0, 0)|0j((^ + $(0, 2k + l)|0)(2lTl| 
+$(2fc+ l,0)|2in)(0| 
+$(2fc-F l,2fc-H l)|2XTl)(2fcTl| , (38) 
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with 



(39) 



and pij :~ {i\p\j). In order to express these formulae in 
the computational (Z) basis, we apply Eq. ([7]) yielding 



D-l 



Pi] 



il—raj 



(40) 



l,m=0 



Therefore, 



^(^' y) = hY. P/™(-l)"'"'™A(; ~m,D), (41) 



where 



n \ sm - 

s{l — ra) 2 



s— — k 



7T{l — ni) 



(42) 



D 



IS 



and Z) = 4fc + 2 as usual. The formula of Eq, 
crucial for our purposes. In fact, it will enable us to 
test the correcting performance of our minimal phase- 
code {|0), |2fc + 1)} in a quantum communication sce- 
nario where the prevalent effect of decoherence is ascrib- 
able to phase damping. 



III. PHASE DAMPING CHANNEL FOR 
QUDITS 

The phase damping (or phase flip) channel for a qubit 
can be defined by the following Kraus decomposition [l^ 



with Kraus operators 



Eo = sl'-±^I,E, = J'-^Z^ 



(43) 



(44) 



where / is the two-dimensional identity operator and Z is 
given by Eq. (U with D = 2. One can describe the phase 
damping channel in an equivalent way, by adopting two 
different Kraus operators, related to those of Eq. (|44l) by 
a unitary transformation 

£{p) ^ E'.pE'^ + E[pE'^ , (45) 

where now 

i?^ = |o)(o| + r?|i)(i| , i?; = vT^|i)(i| . (46) 

The two Kraus decompositions of the phase damping 
channel for qubits, Eq. and Eq. (|^5|) . suggest two 



different generalizations to the general case of dimension 
D. The decomposition of Eq. can be straightfor- 
wardly generalized as 



D-l 



(47) 



m=0 



with 



E„,, — 



D-l \ fl-vY A + 77 



(48) 

which can be seen as a particular example of a Weyl 
channel [13], which is generally defined as 

D-l 

p ^ £{p) = (^"^'") P (^'"^")^ , (49) 

m,n— 

with < n,n,n < 1 such that J2ti^,n=0 ^rn,n = 1- 

There is a however a different way to define the phase 
damping channel for a _D-dimensional spin system, which 
is more closely related to the usual physical meaning of 
phase damping. In fact, phase damping usually means 
that decoherence affects the elements of a given basis 
leaving the elements of the complementary basis un- 
changed. As a consequence, in the basis of the unaf- 
fected states, decoherence destroys only the off-diagonal 
elements of the density matrix of the qudit state. (It 
is evident that here we are conventionally fixing the ele- 
ments \j) of the computational Z-basis as the unchanged 
ones, but it is understood that a complementary damping 
channel can be symmetrically defined). 

The second definition of the phase damping channel for 
a qudit is suggested by Ref. [13| that has shown that the 
Kraus decomposition of Eq. (|45p is equivalent to another 
Kraus decomposition. 



£{p)^J2^^P^^^ 



(50) 



1=0 



with an infinite number of Ei, which in the case of qubits, 
are given by 

. , , , , r7(— 21n77)*/^ , , , , , ^ 

= <5.o |0) (0| + j=l^ |1) (1| . (51) 



The straightforward generalization of such channel in D 
dimensions has been already studied in Ref. [TJ] and has 
Kraus operators 



D-l 



In all these examples, the parameter rj e [0, 1] describes 
the strength of the damping [3]. Such a parameter 
can be assumed independent from D because it usu- 
ally depends only on the bath characteristics, e.g., in the 
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Markov approximation it depends upon the spectral den- 
sity of the bath only [Toj . One can parameterize 77 = e~'^ , 
where 7 is proportional to the probability of a phase er- 
ror, so that phase damping is larger for smaller rj. In 
particular, rj —^ 1~ and 77 — > 0"*" correspond to the weak 
and strong damping limit, respectively. 

It is easy to check from Eqs. ([50]) and (15^ . that an ar- 
bitrary density operator p ^ 
into the output density operator given by 



Ez^iio^'y N) 01 is mapped 



D-l 



(53) 



i.e., we have, as expected, partial suppression of only the 
off-diagonal matrix elements of the state in the compu- 
tational Z-basis. This latter equation shows why this 
second definition of the phase damping channel for qu- 
dits reproduces the usual phase decoherence effect. 



IV. INPUT-OUTPUT FIDELITIES 



The effects of the phase damping Wcyl channel and of 
the conventional phase damping channel on this logical 
state can be described with a unified formalism. In fact, 
one can write the output state of the channel for the two 
cases as 

1 ^"^ 

£ W. <^)] = E "'™/-(^' ^ - ^) 10 , (58) 



i,m=0 



where r — \ refers to the conventional phase damping 
channel, r = 2 to the Weyl channel and 



f2{vJ~ m) 



V 



(59) 



1-rj 



I +11 



D-l 



(60) 



All the results can be expressed in terms of these two 
functions associated to each channel. In order to estimate 
the decoherence effects we compute the fidelity between 
the input and output states 



Once we have defined the two kinds of phase damp- 
ing channel for a qudit, the Weyl channel of Eqs. (HT)) 
and (j48p and the conventional phase damping channel 
of Eq. (j53p , we can consider their action upon our mini- 
mal phase code {|0), \2k + 1)}. Here, we encode a logical 
qubit into a qudit (with dimension D = Ak + 2) hy means 
of these codewords, and we analyze the effects of phase 
damping with and without error recovery. Such effects 
are quantified in terms of fidelity of the output logical 
state with respect to the input. The results are then 
compared to the case where a bare qubit is sent through 
the channel, i.e., when neither encoding nor decoding is 
performed. 

Let us encode an arbitrary pure state cos(6'/2) |0) + 
e"^sin(0/2) |1) of a qubit into a coherent superposition 
of phase-codewords, i.e.. 



\e,(j)) =cos-|0)-f-e*'^sin-|2fc+l) 



(54) 



In the computational basis, the logical state of Eq. ([54]) 
reads 



D-l r 



E 

^ 1=0 



COS- + (-l)'e"^ sin- 
2 ^ ^ 2 



10 , (55) 



and the corresponding density operator is given by 

D-l 

D 



D-l 

p(0,<^) = i0,0)(0,0i = - E 



(56) 



i,m=0 



where 

D.ira := [cos ■ 



ife^'l' sin [cos | + (-l)"e-'* sin |] 



(57) 



D-l 



^ E fr{ll,l-m) 



(61) 



l.m=0 



and, then, average this quantity over all the possible in- 
put states 



^darrip • — . / SlTl 

47r Jo 



dcj) Fdamp{d, 4>) 



D-l 



^ E [3+(-l)'-"]/.(r],^-m) 



(62) 



Lm=0 



The behavior of the averaged fidelity is shown for the two 
cases in Fig. El where (a) refers to the r = 1 conventional 
phase damping channel and (b) to the r = 2 Weyl chan- 
nel. In both cases we note that the decoherence effect of 
the channel increases with the dimension D. 

Let us now apply the recovery map of Eq. ([55)1 to the 
corrupted state of Eq. ((58)) . The output (recovered) state 
is then given by 



PreciO,(f>) ■.^£B,{£[p{0,m 



(63) 



By inserting the matrix elements pim = D ^^imfriVi 
m) of the corrupted state £[p{9, </>)] into Eq. (gT]) we get 
the corresponding $ coefficients 



D-l 



^^^^y) = ^ E ^iMv,l-m){-ir'-y^A{l-m,D) , 



Lm=0 



(64) 

which, in turn, must be substituted into Eq. psp in order 
to give the final explicit expression of the recovered state 
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FIG. 2: Averaged input-output fidelities Fdamp in the case of 
a damped qudit. (a) refers to tfie conventional phase damp- 
ing channel, while (b) refers to the Weyl channel. Fidelities 
are plotted versus the damping parameter 77 and for different 
dimensions D = 30, 18, 6, 2 from bottom to top. 



Prec{Q,4')- The input-output fidelity then becomes 

1 



cos' - + i-iy 



, smf 



+ (-l)'^e-^ + (-l)'-™sin^^ 



(65) 



and its average over all possible input states takes the 
form 

Free ■■= -r / Sin^de / d(b FreeiO,^) 

47r Ja Jo 

D-l 

= ^ E [3+(-l)'-"]/.(»L;-m)Aa-,7i,Z?) , 



Lm=0 



(66) 



Note that the recovery fidelity of Eq. has the same 
form as the damped fidelity of Eq. ([62]) except for the 
presence of the kernel-like term A{1 — m, D) in Eq. 
This term formally takes the recovery operation into ac- 
count and is therefore responsible for the very different 



behavior of Free and Fdamp- Such a term is equal to 1 
only in the trivial case of a qubit (£) = 2 <^ fc = 0) for 
which we have Free = Fdamp-, as is intuitively expected. 
In fact, in this case, the logical qubit is simply encoded 
into another qubit and therefore remains unencoded. 

Our qudit phase-code can be compared with a n qubit 
repetition code. Actually, for a given D = Ak + 2, the in- 
teger n should be chosen as odd [log D~\ , that is as the odd 
integer closest to logi? from above. For a n qubit repe- 
tition code, starting from Eq. (P^ . we straightforwardly 
get 

(n-l)/2 , X /, , X ri-fc /, x fc 

p _ I n\ fl + ri\ fl-V 

^ rec / ^ 



k \ 2 



1 " 



3 ^ V fc / V 2 



(67) 



Notice, that this result holds for both channels, and 
moreover for n = 2 it corresponds to one qubit code, 
i.e. to unencoded qubit. 

The correcting power of our qudit phase-code is shown 
in Fig.[3l where the recovery fidelity of Eq. ((66)l has been 
plotted as a function of the channel decohcrencc param- 
eter rj for several dimensions D = 4fc + 2, for the two 
examples of phase-damped channels, the conventional 
phase damping channel in Fig. [3^ and the Weyl chan- 
nel in Fig. [3Jd. 

In Fig. [3^ we see that in the case of conventional phase 
damping the qudit code performs better and better for 
increasing dimension and always outperforms the unen- 
coded case D = 2 (lower dashed line in Fig. [3^). More 
precisely, in both limits ^ 1^ (weak damping) and 
77 — > 0"*" (strong damping) the correction scheme does 
not depend upon D. However, the improvement with 
respect to the unencoded transmission of the qudit is re- 
markable in the intermediate regime (see Fig. [3^). We 
have also made the comparison with the block-encoding. 
Actually, the dotted lines from top to bottom in Fig. [3^ 
refer to repetition code of dimension 32, 8, 2 (using 5, 3, 1 
qubits) respectively. This code always outperforms the 
qudit code proposed here for any D, showing that even 
though useful, our qudit codes are not optimal in the 
case of the conventional phase damping channel. This 
is however not surprising because the qudit code is very 
different from repetition codes and it is not designed to 
cancel errors at first order in the error probability as do 
the latter codes. 

The situation for the Weyl channel shown in Fig. [5]d 
is more involved. In this case, the correcting power of 
the qudit code increases for increasing D only at small 
phase damping 77 — > 1, while worsening for increasing di- 
mension in the strong damping limit 77 ^ 0. This means 
that for a Weyl channel the qudit code outperforms the 
unencoded case D = 2 only at large rj {rj > 0.7) and 
therefore it is useful only in the weak damping limit. 
In this limit however, contrary to what happens for the 
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conventional phase damping channel, the qudit code be- 
comes particularly useful because it can outperform even 
the repetition code (dotted lines from top to bottom in 
Fig. [313 refer to repetition code of dimension 32,8,2 re- 
spectively). In particular, while the qudit code of D = 6 
does not outperform the 3-qubit repetition code, the qu- 
dit code of D = 18 outperforms the 3-qubit repetition 
code, and quite remarkably the qudit code of Z? = 30 
outperforms the 5-qubit repetition code. 

Finally, the qudit code turns into a non effective code 
while decreasing 77 (worsening for increasing dimension), 
because it is tailored to correct errors of weight up to 
k while in such a limit errors of higher weight becomes 
more and more probable. 
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FIG. 3: Averaged input-output fidelities Fr^c in the case of a 
damped qudit in the presence of the recovery stage, (a) refers 
to the conventional phase damping channel, while (b) refers 
to the Weyl channel. Fidelities are plotted versus the damp- 
ing parameter 77 and for different dimensions D — 30, 18, 6. 
Dashed lines refer, from top to bottom, to repetition code of 
dimension 32,8,2 (using 5,3,1 quibts respectively). In (a) 
the qudit code is always worse than the repetition code, but 
it performs better than the unencoded case except for rj —> 1 
(see inset). In (b) the qudit code is effective at small phase 
damping ry — > 1, where it can outperform the repetition code 
(see inset). 



A. State-dependent fidelity 

Besides the average fidelity, it is also interesting to an- 
alyze the state-dependent fidelity Frec(^, 0) of Eq. 



Both phase damping channels act on the phase of the 
states and therefore in both cases the eigenstates of Z are 
unaffected by decoherence, as can be easily checked. It 
is interesting to see what happens to the encoded states 
in the presence of the recovery procedure of Sec. IID. 
By construction, the two quantum codewords |+) := |0) 

and |— ) := \2k + 1) of our phase code are stabilized by 
and connected by the logical phase gate Z := Z"^^^^ . 
As a consequence they are significantly affected by both 
phase damping channels even in the presence of error cor- 
rection. For the conventional phase damping channel the 
two fidelities satisfy the simple relation 



F,ec(0, 0) = F,ec(^, 0) ^ i for 7/ ^ 



(68) 



i.e., they are completely dephased under the effect of 
strong phase damping. In the Weyl channel case the ef- 
fect is even stronger and for strong phase damping 77 — > 
one has Frec(0, 0) = Frecij^^^) ~ at large enough D. 
However, one can still find an encoded basis which is un- 
affected by phase damping in the presence of error cor- 
rection. Such a basis is formed by the rotated codewords 
given by the two simple superpositions of the initial code- 
words, 



ICo) 



ICi 



!+) + !-) 

V2 



2k 



(2fc + l)-i/^5^|277) 



(69) 



n=0 
2k 



^-=^ = (2fc + l)-V2^ 127^ + 1) (70) 



One can easily check that these are the eigenstates of the 
logical phase gate Z := Z^fc+i g^j.g connected by the 
single shift operator X, i.e.. 



and 



z ICo) = ICo) , z|Ci) = -|Ci) , 



X ICo) = ICi) , X ICi) - ICo) 



(71) 



(72) 



In other words, these rotated codewords {|Co) 1 iCi)} be- 
have like the Z-cigenstatcs |0) and |1) of a qubit. In 
fact, the corresponding fidelities satisfy, for both phase 
damping channels. 



Free (±7r/2, 0) 



1 



V/c, r] 



(73) 



which is exactly the behavior of the Z-basis eigenstates 
|0) , |1) under the action of phase damping (i.e., they re- 
main unchanged). However, let us remark that these 
correspondences with the single qubit eigenstates hold 
if and only if the qudit codewords are subject to error 
correction. 

One can explicitly show that the two codewords 
{ICo) I ICi)} perfectly restored by error recovery for 
any value of 77 by making use of Eq. (|38|) and Eq. (|41|) . 
from which one can verify that 

$(0, 0) = (-1)"$(0, 2k + 1) = (-1)"$(2A: 4- 1, 0) 

= l>(2fc-M,2fcH- 1) = 1/2 , (74) 
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which imphcs 

iRmcu) iCui)] = \Cu) {Cu\ ■ (75) 

As a final remark, note how the e rror correcting prop- 
erties of the orthogonal choices |+), |— ) and |^o) , ICi) are 
perfectly the same in the regime of weak damping (where 
they correct up to k phase-shifts with exactly the same 
ability), while their performances dramatically split when 
the phase damping becomes heavier and it is no-more re- 
ducible to the standard QEC regime. 

V. CONCLUSION 

In conclusion, we have addressed the problem of how to 
profitably exploit the extra space available by embedding 
a quantum system into a "larger" one (qudit-encoding). 
Such an approach can be useful from the point of view of 
the experimental feasibility of quantum error correction 
schemes, since the dimension D of the encoding Hilbert 
space remains reasonably low. In particular, we have 
considered the minimal D which enables the construc- 
tion of qudit-codes able to restore a logical qubit in spe- 
cific decoherence models. These minimal codes are then 
proven to be efficient in protecting quantum information 
against the detrimental effects of phase damping. This 
study could shed further light into the role that Hilbert 
space dimensions play in quantum error correction. The 



opposite problem of quantum data compression could be 
considered in the same light. That is, data compression 
from Hf" to Hd with 2" > 13, as it has been considered 
in Refs. [l^ and [2l|. 

Possible experimental implementations of these codes 
and the corresponding recovery operations require the 
ability to efficiently implement the generalized Pauli op- 
erators X and Z in an effective D-dimensional system. 
An interesting opportunity is provided by ring-shaped 
optical lattices, which have been proposed as a possible 
quantum simulator of periodic one-dimensional quantum 
systems [l^]. If we place a single atom in a ring-shaped 
lattice with D sites, the ground states in each site are the 
basis states. As a consequence, the amplitude shift X is 
realized by tunneling, while the phase shift Z could be 
realized by applying controlled local Stark shifts to the 
atom. 
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